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Abstract. For many years, there seemed to be significant differences between the contin- 
uum intensity distributions derived from observations and simulations of the solar photo- 
sphere. In order to settle the discussion on these apparent discrepancies, we present a de- 
tailed comparison between simulations and seeing-free observations that takes into account 
the crucial influence of instrumental image degradation. We use a set of images of quiet Sun 
granulation taken in the blue, green and red continuum bands of the Broadband Filter Imager 
of the Solar Optical Telescope (SOT) onboard Hinode. The images are deconvolved with 
Point Spread Functions (PSF) that account for non-ideal contributions due to instrumental 
stray-light and imperfections. In addition, synthetic intensity images are degraded with the 
corresponding PSFs. The results are compared with respect to spatial power spectra, inten- 
sity histograms, and the centre-to-limb variation of the intensity contrast. The observational 
findings are well matched with corresponding synthetic observables from three-dimensional 
radiation (magneto-)hydrodynamic simulations. We conclude that the intensity contrast of 
the solar continuum intensity is higher than usually derived from ground-based observations 
and is well reproduced by modem numerical simulations. Properly accounting for image 
degradation effects is of crucial importance for comparisons between observations and nu- 
merical models. It finally settles the traditionally perceived conflict between observations 
and simulations. 
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1. Introduction 

The contrast of the continuum intensity orig- 
inating from the low solar photosphere rep- 
resents one of many diagnostic tests that 
are needed to check the realism of numer- 
ical models of the solar atmosphere, (see, 
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e.g., iDeubner & Mattijd 1 1975b . In the past, 
there seemed to be a disturbing discrepancy 
in granulation contrast measurements derived 
from observations of the quiet Sun and from 
three-dimension al numerical models (see, e.g., 
lNordiundl ll984). The historic evolution of the 
co ntrast measuremen ts is illustrated in Fig. 1 
by iKiselmanl (2008) for a continuum wave- 
length of 500 nm. Simulations produce values 
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a: orginal observation 



b: deconvolved observation 
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d: original simulation 
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Fig. 1. Comparison of exemplary intensity maps: a) original observation (close-up), b) deconvolved obser- 
vation, c) degraded simulation, d) original simulation. 



between 18 % and 27 %, whereas uncorrected 
observations have contrast s of only a few per- 
cen t points. See also, e.g.. ISpruit et al] d 19901) 
and lSanchez Cuberes et al. I (120001) . 



This apparent discrepancy is now proven 
to be caused by the problematic correction of 
the observed intensity images. Observations 
have to be corrected for degradation due to 
instrumental effects and - in case of ground- 
based observations - the terrestrial atmosphere 
("seeing"). Unfortunately, the properties of the 
degradation are difficult to be determined and 
therefore often only poorly known. The result 
is an often incomplete correction and a corre- 



spondingly too low empirical granulation con- 
trast. 



This short article summarises some 
of the main points from a recent more 
extensive study that we published in 
Wed emever-Bohm & Rouppe van der Voortl 



(120091) . iy using space-borne observations 
with the So l ar Op ti cal Telescope (SOT , 
i Tsuneta et ail 120081; llchimoto etail [20081; 
Isuematsu et ail l2008t IShimizu et al.1 (20081) 
onboa rd the Hinode spacecraft dKosugi et al.1 
2007), only the instrumental degradation 
had to be corrected for, which enabled a 
meaningful comparison of observations and 
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hydrodynamic 
ried out by 



simu l ation, which was car- 
Steffenj d2007l) with C0 5 BOLD 
20081) . For comparison, addi- 



dFrevtag et all 

tional si mulations by Schaffenberge r et aD 



Fig. 2. Intensity distribution of the exemplary in- 
tensity images: original observation (thin dashed), 
deconvolved observation (thin solid), degraded sim- 
ulation (thick solid), original simulation (thick dot- 
dashed). 



state-of-the-art numerical simulations of the 
solar photosphere. 

2. Observations 

We selected 584 images taken with the 
Broadband Filter Imager (BFI) of the Solar 
Optical Telescope (SOT). The blue (450 nm), 
green (555 nm), and red (668 nm) wavelength 
bands have been considered. Only quiet Sun 
regions at positions from disc-centre to the 
limb were chosen. 

SOT has a spatial resolution of ~ 0'.'2 - 
0'.'3, which is clearly sufficient to resolve the 
solar granulation pattern. In a precursory study, 
statistically representative point spread func- 
tions (PSFs) were construc ted for each of the 
chann els considered here (Wedem,eyer-B6hm 
2008). For those, a total of 70 images of 
eclipse and Mercury transit observations were 
employed. The PSFs included the effect of the 
SOT aperture and an approximation for non- 
ideal stray-light contributions. 



3. Synthetic intensity 

Most of our analysis was based on three 
snapshots taken from a recent 3D radiation 



(2006), Wed emever et all (120041) . and 
IStein & Nordlundl d 1998b were used. The 
individual models differ in horizontal grid 
spacing and extent, the frequency-dependence 
of radiative transfer, and the inclusion of 
magnetic fields. 

The intensity synthesis code Linfor3D (see 
http://www.aip.de/~mst/linfor3D_main.html) 
was used to calculate synthetic continuum 
intensity maps for each wavelength band 
and all simulation snapshots for heliocentric 
positions from disc-centre to limb with an 
increment of A/i = 0.05. 

4. Comparison of observed and 
synthetic images 

We now compare exemplary intensity maps 
for the blue channel. Next a close-up from 
the original filtergram, we show the corrected 
intensity image in Fig. [T] The image degra- 
dation due to the instrument was corrected 
for by deconvolution with the detailed PSF 
In this example, the granulation contrast in- 
creases from 12.4% to 25.4%. This compares 
to 25.0% for the original simulation snap- 
shot. Degrading the latter by applying the PSF 
decreases the synthetic contrast to 12.9 %. A 
more detailed mean of comparison are the in- 
tensity histograms displayed in Fig. [2] The 
original simulation has a much broader dis- 
tribution than the original observation and ex- 
hibits two distinct peaks at darker and brighter 
than average intensity values, which is not dis- 
cernible from the narrow observed distribution. 
Correcting for the influence of instrumental 
image degradation, however, produces a close 
match of observed and synthetic intensity dis- 
tribution. 

5. Conclusions 

We conclude that the traditionally perceived 
conflict between observations and simulations 
in terms of granulation contrast can be dis- 
missed. Modern simulations are obviously suf- 
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ficiently realistic to reproduce the main char- 
acteristics of the (lower) solar photosphere in 
quiet Sun regions. It is however important to 
note that both sides have small but noticeable 
intrinsic variations, e.g., caused by selection of 
the field of view. One can therefore not expect 
an exact match of the contrast numbers. 

The granulation contrast as a single num- 
ber is a poor test of the realism of simula- 
tions. Evaluating a complex phenomenon like 
the solar surface convection requires more sub- 
stantial tests. In this respect it is essential to 
note that also the centre-to-limb variation of 
the contrast and the power spectral density of 
the continuum intensity from observations and 
simulations are in good agreement. 

Detailed PSFs with reliable estimates 
of the stray-light contributions are crucial 
for quantitative comparisons as shown here. 
Ideally one would measure the PSF exactly 
simultaneous to recording an intensity im- 
age. As this is essentially not possible, a 
robust correction for the influence of image 
degradation requires significant statistics 
employing a l arge number of inte n sity i mages 
as done b y Wedemever-Bohm (200 8J) a nd 
Wedemeyer-Bohm & Rouppe van der Voort 
(2009). Studies using PSFs and intensity maps 
based on a single or a few images only are 
potentially misleading. 

The next step towards realistic models con- 
cerns the refinement of the thermal structure 
and velocity field above the lower photosphere. 
An obvious way to test it are detailed com- 
parisons of spectral lines as they are cur- 
rently employed in the c ontext of abundance 
determinations (se e, e.g.. lAsplund et alj|2000b 
ICaffau et alJl2008h . 
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